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В существующую аналитическую модель петли гистерезиса, заданную в параметрическом виде, 
внесен ряд улучшений. В частности, в модель введены три фазовых сдвига, которые позволяют 
плавно наклонять петлю гистерезиса в точке расщепления на требуемый угол, а также плавно 
изменять кривизну петли. В результате погрешность аппроксимации петли гистерезиса улучшен- 
ной моделью не превышает 1%, что в несколько раз меньше погрешности аппроксимации петли 
существующей моделью. Улучшенная модель способна аппроксимировать большинство из из- 
вестных типов статических симметричных петель гистерезиса, встречающихся в практике физи- 
ческих измерений. Модель позволяет строить гладкие, кусочно-линейные, гибридные, частные, 
зеркально-отражённые, обратные, реверсивные, двойные и тройные петли. Одним из возможных 
применений разработанной модели является линеаризации пьезосканера зондового микроскопа. 
Улучшенная модель может оказаться полезной при решении задач имитационного моделирова- 


ния научных приборов, включающих в себя элементы с гистерезисом. 
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|. ВВЕДЕНИЕ 

Явление гистерезиса широко распространено в природе, оно часто встречается во 
многих областях науки и техники, включая приборы, используемые в научных исследова- 
ниях. " ^ 3 Существует ряд довольно сложных аналитических моделей, описывающих ука- 
занное явление.” Одной из простых моделей является аналитическая модель, предло- 
женная в работе 1. В данной модели семейство петель гистерезиса описывается сле- 


дующими параметрическими уравнениями 


х(&) = асоз" и $” а, 


(1) 


у(а) =6 эта, 


где а - действительный параметр (0=0...2х); а — х-координата точки расщепления (см. 
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у(2) 
1 = 


|- Рис. 1); б», 6, —- координаты точки насыщения; т — поло- 
жительное целое нечеётное число (т=1, 3, 5, ...), опреде- 
т ляющее кривизну петли гистерезиса; п — положительное 
целое число, определяющее тип петли гистерезиса и её 
кривизну. При пП=1 возникает петля типа Лист, при 
в П=2, 4, 6,... — петля типа Месяц (Бумеранг), при 


П=3, 5, /,... - петля типа Классическая. При возрастании 











— РЬ 


параметра а движение по петле происходит в направле- 


| | | | | 
0.5 0 0.5 1 х( <) 


Рис. 1. Петли гистерезиса типа нии против часовой стрелки, при убывании — по часовой 
Лист (пП=1), Месяц (Бумеранг, 

п=2) и Классическая (п-3). Пло- СТРелке. Начало петли (@=0) и конец петли (а=2л) нахо- 
щадь всех трех петель одинако- 
вая. 





дятся в точке расщепления а. 

Главной отличительной чертой модели (1) является 
её простота. Модель интуитивно понятна, она позволяет быстро создавать петли гисте- 
резиса требуемого типа и легко определять ключевые параметры этих петель а, 6,, 6, т, 
п. Для большинства практических задач точности аппроксимации модели (1) 1.5-6%' 
вполне достаточно. Однако встречаются случаи, когда требуется более высокая точ- 


> способна аппроксимировать петли 


ность. Предлагаемая в статье улучшенная модель 
гистерезиса с погрешностью 1% и менее. 

Модель (1) охватывает большую часть известных типов статических симметричных 
гладких петель гистерезиса. Улучшенная модель позволяет с более высокой точностью 
контролировать наклон и кривизну гладких петель (см. секцию П.А.2). Помимо гладких пе- 
тель улучшенная модель позволяет строить различные кусочно-линейные петли (см. сек- 
цию [1.В), а также гибридные петли (см. секцию 1|.В.1.0), в которых прямолинейные участ- 
ки сочетаются с криволинейными. Кроме того, с помощью улучшенной модели из гладких, 
кусочно-линейных и гибридных петель, а также из их сочетаний можно создавать прочер- 
чиваемые непрерывно двойные (см. секцию 1.С) и тройные петли (см. секцию 1.0), как 
самопересекающиеся так и самонепересекающиеся. В разделе |.Е даны формулы, по 
которым можно определять площади петель гистерезиса. В целом в статье представлен 
общий подход, пригодный для аппроксимации большого числа самых разнообразных пе- 
тель гистерезиса. 

Для облегчения построения, анализа и идентификации рассматриваемых петель гис- 
терезиса к статье прилагается дополнительный материал в виде рабочих листов Матка- 
да® (МайФонН, ЦЗА).5 Чтобы избежать загромождения статьи, в ней не приводится вывод 
некоторых формул. Подробный вывод этих формул можно найти в дополнительном ма- 


териале. 


||. ОПИСАНИЕ УЛУЧШЕННОЙ МОДЕЛИ 
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А. Гладкие петли гистерезиса 
1. Дополнительные представления петли гистерезиса 
а. Представление в виде суммы нерасщеплённой петли и кривой расщепления 
Петлю гистерезиса (1) всегда можно представить в виде суммы двух параметрических 
кривых 
Ха) = х, (а) + х. (а) 
у(а) = у, (в) + у. (а) 


где х!(а=бузт”а, у (а=Б,узта — нерасщеплённая петля; х›(а=асоз$”а, у2(0)=0 — кривая 


(2) 


расщепления. Данное представление удобно использовать при рассмотрении преобразо- 
ваний, изменяющих наклон и/или кривизну петли гистерезиса (1): вначале наклоне- 
нию/искривлению подвергается нерасщеплённая петля, после чего полученный результат 
расщепляется простым добавлением кривой расщепления.° Кроме этого, изме- 
няя/добавляя направление действия кривой расщепления можно создавать двойные (см. 
раздел ||.С.3) и тройные петли гистерезиса (см. раздел П.О.3). 
р. Представление в виде частотного спектра 

Используя формулу Муавра, порождающую функцию х(а) в модели (1) можно также 
представить в виде суммы косинусов и синусов с кратными частотами (частотного спек- 


тра) 


5 
й 


т-1 








а. т: п-1 
х(&) = — УС" соз((т-2К)а) + 5» -У (1) 2 "Ск зт(п-2К)а), 3 
ко 27 ко (3) 
у(а) = 6, эта, 
где С’ — биномиальный коэффициент (К, [ -— целые положительные числа); 


СК = ИК -К)|, если 0<К<|, иначе С*=0 (ниже в (6), (25) С*=0, если К-— действительное 
число). Уравнения (3) справедливы для нечётных п, уравнения для чётных п приведены в 
дополнительном материале. В соответствие с (3) при, например, т=п=З уравнения, опи- 


сывающие петлю гистерезиса Классическая (см. Рис. 1), выглядят следующим образом 


х(е) = [Зсоза+ соз(3а/ |+ т Ззта- з(3а]| 
4 4 (4) 
у(а) =Б ‚эта 
Легко видеть, что запись порождающей функции х(а) в виде (3) фактически является 
разложением х(а) в ряд Фурье по нечётным гармоникам 
| 
х(<) = БА, +У (А, соз(Ка) + В, зт(Ка!) 
К=1 
у(а) =Б эта, 


где коэффициенты Фурье Ах, Бк определяются по алгебраическим формулам 
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т-кК 
Я ЕН 
С... 


т-1 





А, = 


%) 


ь (6) 
к Ь, "и 
в, =). 215% 2. 


Амплитуда постоянной составляющей А (К=0) и амплитуды Аки В+ всех чётных гармоник 





(К=2, 4, 6, ...) в (5) равны нулю. Величина / задается равной наибольшей из степеней ти 
п. Функция Поог в выражении для Вхкв (6) — необязательна, она используется только для 
того, чтобы избежать появления комплексных чисел при чётных К. В качестве примера, 
коэффициенты Фурье порождающей функции х(а) петли Классическая (т=п=3), изобра- 
женной на Рис. 1, выглядят следующим образом: Ах=(0, За/4, 0, а/4), В„=(0, 364, 0, -644). 
Располагая коэффициентами Фурье Ак, Вх (6), порождающую функцию х(а) можно 
также представить в виде 
| 
х(<:) = >. Ат, соз(Ка-— ф, } (7) 


где амплитуды Ать и фазы 9 гармоник определяются по формулам 





Например, петля гистерезиса Классическая (т=п=3), показанная на Рис. 1, имеет ампли- 


туды Ать=(0, З./а“ + 65° 14, 0, \/а’ +6: /4) и фазы 1апф»=(0, Бра, 0, -Буа). 
Порождающую функцию х(а) также можно представить в экспоненциальном виде 
/ : 
Ма» Се *. (9) 
К=-/ 
где /— мнимая единица, а комплексный коэффициент Фурье Ск определяется следующим 
образом 
5(А В, } К=1,3,5,...., 


С, = (А +В.) к=-1-3--5,. (10) 


О, К =0,-+ 2,4, 6,.... 


Располагая комплексным коэффициентом Фурье С», порождающую функцию х(а) так- 


же можно представить в виде 
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/ / 
х(а) = >`Ате"“*® = У` Ат, соз(ка- ф, ) (11) 
К=-[ К=-[ 


где амплитуды Ать и фазы 9 гармоник определяются по формулам 


Ат, =_/Ве?(С, )+ т"? (С, ), 

ф, =- Аго(С, = _ аголат (С) | (=) 
Ве(С, ) 

где Ве, |т и Ага — реальная часть, мнимая часть и аргумент комплексного числа, соответ- 

ственно. 

Представление порождающей функции х(а) в виде частотного спектра (3), (5), (7) и 
(11) позволяет синтезировать петли гистерезиса требуемой формы, наклона и кривизны 
путём изменения амплитуд и фаз гармонических составляющих и путем добавле- 
ния/исключения гармонических составляющих с определенными значениями амплитуд и 


фаз.° Данный подход делает возможным построение гладких петель гистерезиса практи- 





чески любой формы. 


у(<) 
1 | 


На Рис. 2 показан пример синтеза Классической петли 





> д, =1.0 
гистерезиса сложной формы’ из гармонических состав- ь,=10 | 
0.5 ая 
ляющих. Порождающая функция петли х(а) образована е = 03 е, Ата = 0.70 





Атз = 0.23 
Ать = 0.10 - 
Атт = 0.02 

$1 = 73° 

фз = -92° 

ф5 = 98° и 
ф7 = -2° 


путем суммирования первых четырех нечётных гармоник ° 


(т=7). Петля гистерезиса проведена точно через 14 


—0.5 Г 


предварительно заданных точек, шесть из которых - это 


ключевые точки петли +а, +а,, +6 (где а, — вертикальное с. 








| | | | | 
0 0.5 1х( <) 


— -0.5 
Рис. 2. Петля гистерезиса слож- 
фазы гармонических составляющих определены числен- НОЙ формы типа Классическая, 

. _в ПорРождающая функция х(а) ко- 
но в ходе решения системы нелинейных уравнений. торой образована суммирова- 
нием первых четырех нечётных 
гармоник. Петля проведена точ- 
первой гармоники, остальные гармоники на площадь пет- но через 14 предварительно за- 


данных точек. Площадь петли 
ли влияние не оказывают (см. раздел П.Е.1). зависит от амплитуды и фазы 





расщепление, см. разделы 1.С.З, 1.0.3). Амплитуды и 


Площадь петли зависит только от амплитуды и фазы 


2.Использование фазовых сдвигов только первои гармоники. 
Одним из основных изменений, внесённых в ранее предложенную модель петли гис- 


терезиса (1), является введение фазовых сдвигов Ася, Ао, Аз 


х(&) = асоз" (+ да.) + 6, зт"(а- да), 13 
(а) = Ь, зп(а + Аа, ), Сы 


^ 


где а, Б 


Хх 


— скорректированные параметры для а, 6», соответственно. Рассмотрим внача- 


ле действие каждого из трех фазовых сдвигов Ася, Ло, Ао в отдельности. 
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у(а) 
1 


- а. Наклон петли с помощью фазового сдвига Ас 


Фазовый сдвиг Аа. позволяет плавно наклонять петлю 


0.5 


гистерезиса, изменяя угол наклона р=л/2-6 касательной к 
петле в точке расщепления а (см. Рис. 3). В ранее пред- 
ложенной модели наклон петли осуществлялся путём по- 
- ворота системы координат на угол ви предыскажения па- 


—0.5 


раметров петли а, 5», б, путем их поворота в противопо- 


— 











| | | | 





ЛОЖНОМ направлении. ' Более точные формулы наклоне- 


| 
0.5 Л х( а) 


1 0 
Рис. 3. Наклонение петли гисте- вия петли поворотом, позволяющие исключить до этого 
резиса типа Классическая с по- 


МОЩЬЮ фазового сдвига Ла. Имевшее место небольшое смещение точки расщепления 
Наклонение в точке расщепле- 
ния а на угол +9 обеспечивается 
сдвигом на +Лод. Наклонение (Л=А%=А=0)° 
петли вызывает рост её площа- 

ди. Площадь петли с положи- 

тельным наклоном равна пло- 

щади петли с отрицательным 

наклоном. 


— —0.5 


от заданного положения, выглядят следующим образом 


х(а) = х(®)+зтб(Ь, зт0+ Ь, со$ 9)6та— т” а} й 
У(<) = у(@®) + зтб(Ь, соз0- Б, зт@ та - эт” зы 


Из преобразований (14) видно, что наклонённая петля представляет собой результат 
сложения исходной петли с некоторой кривой, обеспечивающей наклонение исходной 
петли (см. раздел 1|.А.4). 

Поскольку внесение фазового сдвига Ао’ ведёт к изменению координат точек расще- 


пления а и насыщения Б,, то требуется коррекция координат этих точек. В улучшенной 


^ 


модели (13) скорректированные параметры а, Б 


, находятся из следующей простой сис- 
темы уравнений, составленной для точки расщепления о=0 и точки насыщения а=л/2 
(Аа=А0%=0) 

Асоз” (0+ ла.) +6. зт"0=а, 
Л (15) 


^ у[й Е: 
асоз"| —+да, |+6, эт” — =Б›, 

2 2 
откуда легко определить 


а 


а= 
0$” Да. (16) 


г _ т 

6, =б +айап” Аа. 
При 7=1 петли типа Лист (п=1) не зависят от фазового сдвига Аа. Фазовый сдвиг Ас 
требуемый для наклонения петли в точке расщепления а (а=0) на заданный угол 9 рас- 


считывается по формуле? 
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(17) 


Нерасщепленные петли (а=0) нельзя наклонить с помощью фазового сдвига Даа. 


р. Изменение кривизны петли с помощью фазового сдвига Лесь 


Фазовый сдвиг До по- 
зволяет изменять кривиз- 
ну петли гистерезиса (см. 
Рис. 4). В отличие от па- 
раметра т фазовый сдвиг 
А > обеспечивает плавное 
изменение кривизны пет- 
ли. Также как и в случае 
фазового сдвига Ао, при 
фазовом сдвиге Ас необ- 
ходимо определить скор- 


ректированные координа- 


^ 


ты а, Б 


х 


точек расщепле- 





-0:5: — 








== 


| | | 
—1 —0.5 0 0.5 








1 х(@) 
Рис. 4. Плавное изменение кри- 
визны петли гистерезиса с по- 
мощью фазового сдвига Ле. 
Площадь петли возрастает при 
увеличении фазового сдвига 


Аср. 


0.5 - 


—0.5 Г 











Я - 


| | | 
1 —0.5 0 0.5 





- 1 х( а) 
Рис. 5. Плавное изменение кри- 
визны петли гистерезиса с по- 
мощью фазового сдвига А. 
Площадь петли уменьшается 
при увеличении фазового сдви- 
га Ао%. 


ния и насыщения, соответственно. Для этого следует решить систему уравнений, состав- 
ленную для точки расщепления а=0 и точки насыщения о=л/2 (А=Ао8=0) 


Асоз" 0+6, зт"(0+ла,)=а, 

‚ ^ 18 

асо$" + Б эт" И Аа, =6, (18) 
2 2 

откуда 


а=а-6 4ап” Да,, 
р 6 


х 


Хх 


^ со" у 
При т=1 петли типа Лист (п=1) не зависят от фазового сдвига До. 
с. Изменение кривизны петли с помощью фазового сдвига Лаз 
Фазовый сдвиг Дов, также как и Ас, позволяет плавно изменять кривизну петли гисте- 
резиса (см. Рис. 5). Составляя уравнения для точки расщепления а=0 и точки насыщения 


о=т/2 аналогично тому, как это делалось выше, получаем (Аби=А=0) 
асоз” (0-—ла.)-+ 6, зт"(0-ла.)=а 


асов" — ле +Ь, ыи" — ле =6, 
2 2 


откуда находим скорректированные параметры 
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асоз” Да, +6 3т” Аб 


а = , 
$1” Аа + С0$""" Аб 


В = ый Аа. — м АО. | 
п” Да. +с0$ ° Да. 
При т=1 форма петель типа Лист (л=1) не зависит от фазового сдвига Дов. 

Действие фазового сдвига Лоз=Аса противоположно действию фазового сдвига А%=Аа 
(сравни Рис. 5 с Рис. 4). При этом петля, построенная при А=Аа, немного отличается от 
петли, построенной при Аоз=-Аа. Различие возрастает по мере увеличения |А^4. Фазовый 
сдвиг Лаз вызывает смещение параметра @ на величину Лоз. В результате в точке рас- 
щепления а параметр оравен -Ась, а в точке насыщения 6 -— х/2-Аов. 

а. Одновременное использование нескольких фазовых сдвигов 

В общем случае, когда одновременно задействованы все три фазовых сдвига Дау, А 

и Аоз, система уравнений, составленная для точки расщепления а=0 и точки насыщения 


о=т/2, выглядит следующим образом 


Асоз” (0+ да, - Аа, )-+ Б, эт" (0+ла, -да.)=а 


й ^^. 22 
асов" ТА - да, | +5, ы + а. -да, |5. (22) 
Решая систему (22), находим искомые скорректированные параметры 
РИ асоз” (ла, - да, )-6эт"(ла, - да.) 
эт” (да — Аа. ) т” (Аа — да, + соз" (ла — да. )соз" (да, - да. (5) 
Е азт” (Аа — Аа.) +6 ‚со$” (Аа, — Да, )} 


х 


` эт” (да — Аа, эт” (Аа - да, )+со$” (Аа, — да. )соз" (ла, - Аа, 
Подставляя полученные параметры а, Ь, в (13), легко показать, что при т=1 петли типа 
Лист (л=1) не зависят от фазовых сдвигов Да, Ась, а форма этих петель не зависит от 
фазового сдвига Аоз. Формулы (16), (19), (21) являются частными случаями формул (23). 
Хотя действие фазового сдвига Аоз противоположно действию фазового сдвига Ао, 
полной нейтрализации одного сдвига другим при их совместном использовании 
ДА 0р=А08=Аа (Аб=0) не происходит. При задании равных друг другу фазовых сдвигов 
А %=А08=Аа их совместное действие противоположно действию фазового сдвига Ас=Ла. 
Поэтому для наклона петли гистерезиса в точке расщепления вместо использования фа- 
зового сдвига Аои=ЛАа можно задавать А%=Аоз=-Аа. Поскольку действие фазовых сдвигов 
Аи А0%=-А схожи, то при практическом применении модели (13) можно ограничиться 
использованием только двух фазовых сдвигов: До и Ло (см. раздел 11) или Да и Ада. 


Формулы для определения скорректированных параметров а и 6, при использовании 


фазовых сдвигов Ася, А получаются из формул (23) при подстановке Ао%=0; при исполь- 
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зовании фазовых сдвигов Ао, Аоз — при подстановке Ао»=0. 


Заметим, что при Аа=А%=А0=Ао=0, где да — произвольное действительное число, 


^ 


скорректированные параметры а, 6, улучшенной модели (13) вырождаются в параметры 
а, 6, исходной модели (1), соответственно. Таким образом, в рассматриваемом случае 
петли гистерезиса, построенные согласно модели (13), по форме совпадают с петлями 
модели (1), но только их параметр а оказывается сдвинутым на величину да. Вывод о 
вырождении модели (13) в модель (1) также вытекает из следующего рассуждения. По- 
скольку одновременное действие равных друг другу фазовых сдвигов А%=Ааз=Ла проти- 
воположно действию фазового сдвига Ая=Аа, то при одновременном задании трех оди- 
наковых фазовых сдвигов Ао=А%=Аоз=Аа должна возникнуть петля (1). 

Также как и исходную модель (1) улучшенную модель (13) можно представить в виде 


суммы гармонических колебаний 








т-1 п-1 
2. Е п-1 
х(а)=—9 - У Ск с0$((т-2К(а+ Да, )) + >. -У (1) 2 "Ск зт(п-2К(а+ да, 
2" К=о 2 К=0 (24) 


У(а) = 6, зп(а + Да, ) 


Из представления (24) можно также напрямую вычленить коэффициенты Фурье Ах, Ве 











т-к к р пк 
А, = — С. бане 5 -С,? эи(КАа,, 
^ т-К м п-К (25) 
а от. Г Ба 
В, ета „2 эп(КАа, )+ (-1): 2 ре соз(Клаь ) 


по которым легко определяются амплитуды Атх и фазы 9х гармоник порождающей функ- 
ции х(а) (см. формулы (8), (12)). 

Согласно определения (13), для получения на выходе уа) гистерезисного элемента 
гармонического сигнала на его вход нужно подать сигнал вида 


х(@)=асоз" (сд) + Б, зт"(о+л ось) (или в соответствии со спектральным представлением 


(24) сумму гармоник (7)). Требуемый начальный сдвиг фазы Аа гармонического сигнала 


^ 


получается путём задания соответствующих скорректированных параметров а, Б 


Хх? 


ВЫ- 


числяемых по формулам (23). Из сказанного следует, что петля гистерезиса (13) облада- 
ет фильтрующей способностью. 
3.Петли гистерезиса типа Летучая Мышь и Астро 


В работе 1 петлю гистерезиса типа Летучая Мышь (Бабочка) получали, взяв порож- 


дающую функцию у(а) по модулю: У(а) = 65 о. Подобного типа петлю также можно по- 


лучить, если в У(а) вместо модуля возвести зпав чётную степень К=2, 4, ...: 
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Х(<) = х(а), 
У(а) =Б эт“ а ро 


На Рис. ба изображена петля гистерезиса типа Летучая Мышь при К=2. При нечетных 





У( @) 
= 1 => 


У(а) Я( <) 
| = ЧЕ 







„.277-° 


с со 
И 
о 
№ 
о о 


< х 


ии 


0.5 — — 0.5 Е — 0.5 —- 


пи 
© 


а - 
— 
сл 
о 


а2 = 0.1 

Ь1к = 0.27 | Ь>х = 0.73 
62, = 0.5 

т5 = 3 

По = 3 





А; = Ар = В. = Вь= 1.0 











аа” К1 = К2= 11 = 12 =1 _ 




















| | | | | | | | | | | | | | | 
—1 0.5 0 0.5 1Х( <) ых 0:5 0 0.5 1Х( <) = -0.5 0 0.5 1х( а) 





(а) (6) Рис. 7. Наклонная Классическая 
Рис. 6. Петля гистерезиса (а) Летучая Мышь (Бабочка), (0) Астро. петля (зелёная) — результат 
степенях К=3, 5, ... получаются петли типа Астро.'3 На СЛОоЖения двух петель гистере- 


зиса — петли 1 Лист (красная) с 


Рис. 66 изображена петля гистерезиса типа Астро при петлей 2 Классическая (синяя). 
Для сравнения пунктирной ли- 
К=3. нией показана петля Классиче- 


4. Арифметические действия над петлями ская, наклоненная в точке рас- 
щепления а на такои же угол 


гистерезиса 6=15° с помощью фазового 
вига Ао. 

Петли гистерезиса (1)-(5), (7), (9), (11), (13) и их ком- бАВИГа ли 

поненты можно складывать, вычитать, умножать на число, возводить в степень, а также 


использовать в качестве аргументов некоторых функций 
— К; 
х(а) — У Ах, (©) 
/ 


27 
уе) = У ву (а) и 


где А, В, — действительные коэффициенты; К, /-— положительные целочисленные степе- 
ни, обычно нечётные. Определенные над петлями арифметические операции (27) позво- 
ляют изменять форму, наклон и кривизну петель гистерезиса, а также строить двойные 
петли (см. раздел |.С.4). Допустимо складывать петли разных типов. 

На Рис. 7 показан пример сложения двух петель гистерезиса разных типов — петли 
Лист и петли Классическая (а›=а-а!, б>.=б»-6ч», б2,у=б,-61,). Параметры результирующей 
петли а, 6,, 6, т, 9 даны в левом верхнем углу рисунка. Следует обратить внимание на 
то, что посредством сложения с петлей типа Лист Классические петли можно наклонять в 
точке расщепления а на угол 9 (р=л/2-0 — угол наклона касательной в точке а). Параметр 
петли 61» (или бу, б>у, б,) находится с учётом заданного угла 9 по формуле? 

6, = (5, + 6, ,\апб = 6, 1апб. (28) 


В общем случае любую петлю гистерезиса можно разложить на бесконечное множество 
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Рис. 8 Наклонение петли гисте- (а) (6) 


| 
0.5 1х( а} —1 -0.5 0 0.5 1х(а) —1 —0.5 0 0.5 1х(а) 


резиса Классическая путём пе- Рис. 9. Изменение кривизны петли гистерезиса Классическая пу- 
рекоса системы координат на тём перекоса системы координат на угол кв направлении оси у. 
угол дв направлении оси х. Пло- Наклонение петли в точке расщепления на угол (а) в=15°, (6) 
щадь всех петель одна и та же 6=-15°. Петли, наклонённые на любой угол 09, у которых все ос- 


при любых углах перекоса 6. тальные параметры одни и те же, имеют одинаковые площади. 


петель, поскольку указан- 
ным выше способом каждую петлю из пары складываемых петель всегда можно предста- 
вить в виде некоторой пары каких-то других петель. 

В дополнительном материале показано, что действие следующего преобразования 
(Аа =Д 0=А 0%8=0) 


х(<) = х(а)+ р, ‘ап 9 И ш” а 


(29) 


У(а) = у(а) 
построенного на перекосе системы координат на угол дв направлении оси х, эквивалент- 
но наклонению петли посредством сложения с петлей Лист (см. Рис. 7). На Рис. 8 изо- 
бражены Классические петли гистерезиса, наклонённые в точке расщепления согласно 
преобразованию (29). 


Дополнительный перекос системы координат на угол кв направлении оси У 


х(<) = х(а) + пап в(Ь. (ап, {па — $1” а) 

У(а)= у(@)+Б. пап «па зт" а р 
позволяет плавно изменять кривизну петли. На Рис. 9 показаны построенные с помощью 
преобразований (30) наклонные Классические петли с разными кривизнами. Преобразо- 


вания (29) являются частным случаем преобразований (30). 


В. Кусочно-линейные и гибридные петли гистерезиса 

Простейший способ построения кусочно-линейных петель в рамках рассматриваемой 
модели заключается в следующем. Период Т=2л изменения параметра с делится на нуж- 
ное количество К интервалов, где К — целое чётное число (К>4). Подставляя в формулу 
(13) вместо непрерывных значений а, значения а, изменяющиеся с шагом Г/К, получаем 


х, у-координаты точек кусочно-линейной петли гистерезиса, которые соединяем отрезка- 
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у(а) 
| 1 


У(а) 


ми прямых. На Рис. 10 дан ”' 
пример кусочно-линейных 
петель трёх типов Лист, “ 
Месяц (Бумеранг) и Клас- 
сическая, построенных 
для К=12 (шестисегмент- +х. 


ная петля, А/2=6). 














Другие способы полу- | | | 


—1 -0.5 0 0.5 1 х(а) | -0.5 0 0.5 1 х( а) 
чения  кусочно-линейных Рис. 10. Кусочно-линейные пет- Рис. 11. Кусочно-линейная пет- 

ли гистерезиса типа Лист (п=1), ля гистерезиса Лист (Люфт без 
петель гистерезиса за Месяц (Бумеранг, п=2) и Клас- Усов, п=1), гибридный Месяц 
сическая (п=3). Число интерва- (гибридный Бумеранг, п=2) и 
лов разбиения К=12 (шестисег- гибридная Классическая (п=3), 
функций синуса и косину- ментная петля). Площади пе- построенные на трапецеидаль- 


тель примерно одинаковые. ных импульсах. Площадь всех 
са в порождающих функ- трёх петель одинаковая. 


циях х(а) и уа) модели (13) на трапецеидальные, тре- 








ключаются В замене 


угольные и прямоугольные импульсы единичной амплитуды и их сочетания. 
1. Петли на трапецеидальных импульсах 

Заменяя синусы и косинусы в порождающих функциях х(а) и (а) модели (1) на трапе- 
цеидальные импульсы единичной амплитуды р. и ре, соответственно, получаем урав- 
нения кусочно-линейных петель гистерезиса, построенные на трапецеидальных импуль- 
сах 

х(а) = ар" а- (Б, —а)тр" а, ы. 

у(а) = 6, пр, <, (31) 
где р‹(а=игр.( ая 1/4); Т —- период импульсов. Выражения, с помощью которых задаются 
трапецеидальные импульсы 1{гр, приведены в работе 1 и в дополнительном материале. 
Наиболее простые петли получаются, когда верхнее А и нижнее 0) основания трапецеи- 
дальных импульсов связаны между собой как О=За (Т=а9+0=4а). В этом случае форма пе- 
тель не зависит от т. Вычитание в (31) расщепления а из 6, позволяет переместить точку 
насыщения 6 из середины горизонтального участка петли, где о= 1/4, в изначально приня- 
тое место её расположения (см. Рис. 11); здесь параметр а принимает значение равное 
ТВ. 

На Рис. 11 изображены петли трёх типов: Лист (Люфт без Усов), Месяц (Бумеранг) и 
Классическая, построенные согласно формул (31) (Р=За). Петля Лист является кусочно- 
линейной. Петли Месяц и Классическая являются гибридными — в них прямолинейные 
участки сочетаются с криволинейными. В случае, когда 6б,=а, петли (31) вырождаются в 


прямоугольную петлю (Неидеальное Реле без Усов). 
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—1 -0.5 0 0.5 1 х(@) —1 -0.5 0 0.5 1 х( <) —1 -0.5 0 0.5 1 х( <) 


(а) (6) (в) 

Рис. 12. Кусочно-линейные петли гистерезиса: (а) Люфт-Реле-Люфт (тетра-линейная петля), (6) 
Люфт-Люфт (три-линейная петля), (в) Люфт-Реле (три-линейная петля), построенные на трапе- 
цеидальных импульсах с использованием фазовых сдвигов. 


а. Кусочно-линейные петли 


С учётом фазовых сдвигов Аа, Ао, Ао уравнения (31) приобретают вид 


х() = атр"(и+ да, +6, тр” (а + да, 2 
у(а)=Ь, пр. (&+ да. } (32) 


^ 


Скорректированные параметры а и Б, определяются с помощью следующих формул° 


РИ а\гр” (с, + Да, — Аа. )-Ь лир" (да —- да.) 
р” (ла — Да, тр” (а, + да, - да, )- пр” (а, + да, - да, гр" (да, - Аа, (33 
, Ь 4гр” (Аа. — Аа. )- ар” (а, + Да - Аа. ) 


* ир”(да, — да, тр" (а, + да, -да,)- пр" (а, + да, — ла, р" (да, - да.) 
где о равно Т/8-Ао-+Аов, Г/В-Аа, 3 ГВ-АаФ+Ась и другим значениям? в зависимости от 


используемых диапазонов фазовых сдвигов Ла, Аа, Аоз. Логика вывода формул (33) 


^ 


скорректированных параметров а и 6, для кусочно-линейных петель (32) такая же, что и 


для случая гладких петель. При %=Т/4 рабочие формулы для определения скорректиро- 


^ 


ванных параметров а и Б, можно получить из формул (16), (19), (21), (23) для гладких 


х 





петель (13) простой заменой функций эт и со$ трапецеи- к. 


1 |=. 





| 
> 


дальными импульсами Тр. игре, соответственно. 


ссо 
< х 


и 


— 


На Рис. 12 показаны кусочно-линейные петли гистере- о.- 





э3 


зиса, построенные с помощью уравнений (32) и скоррек- 
тированных параметров (33) (т=п=1). В качестве примера 
на Рис. 13 представлена петля Люфт-Люфт при различ- 


ных значениях фазового сдвига Ао. Согласно определе- 





ния (32) при подаче на вход кусочно-линейного гистере- +1-, 





| 
1 0.5 0.5 1 х( <) 


= . 0 
ЗИСНогО элемента сигнала ВИДа Рис. 13. Кусочно-линейная петля 


лы ^ И гистерезиса Люфт-Люфт при 
ха)=а!р. (0+Ао)+ 6 1гр. (а+Ааё), на выходе у(а) получа- различных значениях фазового 


сдвига А(%. 
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Улучшенная параметрическая модель для аппроксимации петли гистерезиса 




































У(а) У(а) | У(а) 
у= -10° 
1 — 1 а = 0.4 6 Ь | 1 а = 0.4 _Ь 
Бу = 1.0 бу = 0.8 
Бу, = 1.0 БУ = 1.0 
т = 1 т = 1 
0.5 — 0.5 п=1 — 0.5 п- 1 
Т=8 ре 
Ла1 = - Ла = 0 1 
ое -о0 Ао = —1 -о Да = 1.8 
Лаз = 0 1 
—0.5 —-0.5 —5.5 
4 — ь М Е у Т — 
1 1 —6 = — Да " —6 —“ а д. + Да 
1 Ь 8 2 3 
| | | | | | | | | | | | 
— 0.5 0 0.5 1 хХ(а) — —0.5 0 0.5 1 х(а) — 0.5 0 0.5 1 х(а) 
(а) (б) (в) 


Рис. 14. Кусочно-линейные петли гистерезиса с усилением/ослаблением у (а) Люфт-Реле-Люфт, 
(0) Люфт-Люфт, (в) Люфт-Реле, построенные на трапецеидальных импульсах с использованием 
фазовых сдвигов. 


ются трапецеидальные импульсы с заданным фазовым сдвигом До. 
Чтобы получить кусочно-линейные петли с усилением/ослаблением у добавим к (32) 
дополнительный член (кривую), отвечающий за усиление/ослабление (т=п=1)° 
х(е.) = х(е) 
У(@) = уе) + пап а) в, ир, (а + ла зы 
На Рис. 14 показаны петли с усилением (7>0) и с ослаблением (отрицательным усилени- 
ем 7<0), построенные по уравнениям (34). Петли на Рис. 14 без усиления (7=0) соответст- 
вуют петлям, изображенным на Рис. 12. 
Петли типа Люфт без Усов (Лист) (см. Рис. 11) невозможно наклонить в точке расщеп- 
ления посредством поворота системы координат. Тем не менее, с помощью поворота 


этим петлям можно придать усиление/ослабление у(т=п=1, Д=А%=Ав=0)° 


ха) = а1ап рр, а + (6, — 6, {ап ур. [01 
{апр -—1апу | (35) 
У(а) = ап В(х(а) - ар. а) 
где параметры петли а, 6, 6, Б связаны между собой следующим соотношением: 
{апр=б//(6х-а). Петли гистерезиса, построенные на трапецеидальных импульсах с помо- 
щью уравнений (35), представлены на Рис. 15. 
Применяя разновидность поворота системы координат, найдены следующие преобра- 
зования (т=п=1, А=А=А=0)° 
х(о/) = х(.) 
У(а) = у(а) + аапУ (р, @ р, “-— р. а) (36) 
с помощью которых на трапецеидальных импульсах можно строить петли Люфт-Люфт с 
Усилением без Усов" (см. Рис. 16). 


Ь. Гибридные петли 


Наклонение в точке расщепления посредством фазового сдвига Ася гибридной Клас- 
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Г 2=04 


Бр, = 0.58 
Ь, = 0.79 
рат” 

у= -17° 








0.5 Г 


— = 





` а=04 


Ь, = 0.63 
Ь, = 1.0 
В= 77° - 


о 
и = 
/ 








а = 0.4 


Ь, = 0.4 
БУ = 0.82 
05- Д= 90° 
7 = > 








а = 0.4 

Б» = 0.4 

Бу = 1.0 

ДВ = 90° — 
у = 0° 





МЕ 1 ЕЕ - 




















| | | | | 
—1 —0.5 0 0.5 1х( <) —1 0.5 0 0.5 1Х( —1 —0.5 0 0.5 1х( <) 


(г) (д) (е) 

Рис. 15. Кусочно-линейные петли гистерезиса: (а) Люфт с Усилением без Усов (параллелограмм- 
ная, билинейная, упругопластическая петля), (6) Люфт с Ослаблением без Усов, (в) Люфт без 
Усов, (г) Неидеальное Реле с Усилением без Усов, (д) Неидеальное Реле с Ослаблением без 
Усов, (е) Неидеальное Реле без Усов (прямоугольная петля), построенные на трапецеидальных 
импульсах. 











сической петли (п=3, 5, ...), показанной на Рис. 11, не применимо, поскольку при таком 
наклонении на криволинейных участках петли возникают изломы (скачки производной), а 
порождающие функции х(@), у(а) при одних и тех же значениях параметра а приобретают 
горизонтальные участки. Получить петлю с заданным наклоном Д=л/2-0 в точке расщеп- 
ления, сохранив плавность криволинейных участков и горизонтальность прямолинейных 
участков, позволяет следующее преобразование (1=1, Ася =Д%=Аов=0)° 

х(а) = х(а) + Ь, эт (пр. «— р” а} 

У(а) = у@) 


Преобразование (37) для гибридной петли подобно преобразованию (29) для гладкой 


(37) 


петли. На Рис. 17 показаны гибридные классические петли гистерезиса без усов, накло- 
нённые согласно преобразованию (37).'" "*. "3 

Путем применения ряда линейных преобразований перекосаб получена универсаль- 
ная формула (т=1, Ам=А%=А0=0) 

х(а) = ха) +1апб[, —а)\апк- 6, апу+ Ь, пр. а- пр” а) 


38 
у(а) = у(а) + [6 _а\4апк- 6, {ап Ур. а- пр. &)+ атап Уир. сре и -—1тр” а] Е 


Сорупа!1 © 2020 Р. В. Лапшин. Все права защищены 15 


Улучшенная параметрическая модель для аппроксимации петли гистерезиса 


которая позволяет строить 
на базе трапецеидальных 
импульсов гибридные 
петли без 

Рис. 18), 


заданный на- 


классические 


усов"? 13 


(СМ. 
имеющие: 
клон д в точке расщепле- 


(р=л/2-6), 


наклон 


НИЯ заданный 
(усиле- 
ние/ослабление) у прямо- 


линейного участка, тре- 





| а=0.4 
Ь„ = 0.6 
Бу = 0.87 
В = 77° 
и = 17° 





— | = 
| | | | | 
— 0.5 0 0.5 1Х( 1) 


Рис. 16. Кусочно-линейная пет- 
ля гистерезиса Люфт-Люфт с 
Усилением без Усов, построен- 
ная на трапецеидальных им- 
пульсах. 











буемую кривизну к криволинейного участка. Введение 


помимо п=3, 5,... дополнительных параметров К=2, 4, ... 





а = 0.4 
Ви 07 
Бу = 1.0 
й=Э 
Аал = 0 
Аа = 0 
Лаз = 0 


0.5 — 





—6 
| | | 











—1 -0.5 0 05 15а) 
Рис. 17. Гибридные —Классиче- 
ские петли гистерезиса без усов 
с заданным наклоном р=л/2-0 в 
точке расщепления. Петли по- 
строены на трапецеидальных 
импульсах. Площадь всех пе- 
тель одна и та же при любых 0. 


и х; управляющих кривизной, продиктовано тем, что использование фазовых сдвигов Лось, 


АО приводит к возникновению нежелательных изломов на криволинейных участках пет- 


ли. В формуле (38) вместо 
функции р’ можно ис- 


пользовать р: а, где К- 





любое положительное чис- 
ло. При 7=к=0О формулы 


(38) вырождаются в фор- 


мулы (37). 
При необходимости 
формулы (31)-(38) можно 


выразить через треуголь- 
ные импульсы, используя 


следующее представление 





—. 


@) 


у 


1 | 











з®хзэзг 
пни 
я ам 
о 


—> 
Го 
о 

















(а) 


1 Х(а) —1 -0.5 0 


0.5 1х( <} 


(5) 


Рис. 18. Гибридные Классические петли гистерезиса без усов с 
заданным наклоном 0=л/2-6, усилением/ослаблением уи кривиз- 
ной х: (а) Разные усиления упри фиксированных ри х; (6) разные 
кривизны «при фиксированных ри 7 Петли построены на трапе- 
цеидальных импульсах. 


трапецеидальных импульсов в виде суммы двух треугольных (Д=За, Т=а+0). 


р. (а) = Ш. (&+ а) +. (&— =) 


р. (@) = 1 (&+ =) +. (&— 5) 


2.Петли на треугольных импульсах 


Помимо трапецеидальных импульсов {гр в формулах (13) могут использоваться тре- 


угольные импульсы 11, ''° которые являются частными случаями трапецеидальных (6-0, 


=) 


16 
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х(&) = а" (и + да) + Би" (и + Аа), 


у(а) = Б, 1. (© + Аа,). (40) 


Рабочие формулы для определения скорректированных параметров а и 6, получаются 


из формул (16), (19), (21), (23) для гладких петель (13) простой заменой функций $ и со$ 
треугольными импульсами 11. и Шс, соответственно. 


На Рис. 19а представлена кусочно-линейная петля Лист (т=п=1), построенная по 





формулам (40) при различ- кх. ие) | | 


ан 1 да,=0.17 Ь = 





ных значениях положи- 
тельного фазового сдвига о.- 
АО (А=А0%=0). При 
Аа:=0.8 петля является 
разновидностью петли 


Люфт с Усилением без 














Усов (см. Рис. 206); при 1’ 4 . 





Аои=1.5 петля является (а) _ (6) 

Рис. 19. Кусочно-линейная петля Лист (двухсегментная петля, 
роны Мали Люфт без Усов) при различных значениях фазового сдвига Ас, 
Люфт без Усов (см. Построенная на треугольных импульсах по (а) базовым, (6) об- 

ратным базовым уравнениям. 
Рис. 20е). Как видно из 


Рис. 19а, с помощью фазового сдвига Асд возможно управлять усилением 7. 


Показанные на рисунке углы Ви у определяются по формулам 


А И 
Ь/-а Ь,-2а\1апу а(ш, Да, -1)+6, \, Аа, 
(41) 
у бу 1апр _ Б/1и, Да, 





ап = = 
Ь +а Б,+2аапВ ам, да +1) +6. Аа 


Используя формулы (41), уравнения (40) можно при необходимости выразить через углы 
В и/или у При Б,=а петля (40) вырождается в разновидность петли Неидеального Реле с 
Усилением без Усов (см. Рис. 20жЖ). 

Подобно тому, как это было выполнено выше для петель (35), построенных на трапе- 
цеидальных импульсах, усиление/ослабление упетель (40) можно изменять посредством 


поворота системы координат (т=п=1, Ди=А%=Аь=0)° 


хе) = 6 и. а+| 2211 ша 
{ап р —{апу (42 
а+6, \апу- 6 
У(а) = Ь, Ш, и + В, Па, С, 
{ап б —{апу 


в котором параметры петли а, 6», 6,, В связаны между собой соотношением {апр=б/(Бх-а). 


Сорупд © 2020Р. В. Лапшин. Все права защищены 17 












ппя аппооксимации петп 


и гистерезиса _ 





В=77°® 


(а) 


а = 0.4 —- 1-Е а=0.4 — 
Ь„ = 0.60 Ьх = 1.0 
Бу = 0.87 Ь, = 0.6 









в= 77° 





р = 17° 


у= -17° 



































0.5 1 х( 





0.5 


—0.5 























0.5 




















0.5 




















| 
0.5 


| | | | | | | | | | | 


1 х( а) 0.5 





—0.5 0 1 х( <) 1 х( а) 


(и) 
Рис. 20. Кусочно-линейные петли гистерезиса: (а) Люфт с Усилением, (6) Люфт с Усилением без 
Усов (параллелограммная петля), (в) Люфт с Ослаблением, (г) Люфт с Ослаблением без Усов, (д) 
Люфт, (е) Люфт без Усов, (е) Неидеальное Реле с Усилением, (ж) Неидеальное Реле с Усилением 
без Усов, (3) Неидеальное Реле с Ослаблением, (и) Неидеальное Реле с Ослаблением без Усов, 
(и) Неидеальное Реле (триггер Шмидта), (к) Неидеальное Реле без Усов (прямоугольная петля), 


построенные на треугольных импульсах. 
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С помощью уравнений (42) можно строить кусочно-линейные петли, показанные на 
Рис. 206, г, е, ж, и, к. 
Другой способ описания петли типа Лист строится на основе функции обратной к (40) 
(т=п=1)8 
х(<) = 6/1. (а + да. ) 
у(@®)=Ь, 5 Ма де, + 1. (&-+ да) | и 


Хх 


Скорректированные параметры а и Ь, в (43) находятся по формулам° 


Би. (а. + Аа — Ао.) 


а а, 
1. (©, + Аа, — Аа, + Ш. (<, + Аа, -Аа,) 
(44) 
м Би, (а, + Да, - Ао) 
У чи. (а. + да. — да. м, (Ао — Аа) + Ш. (<, + Аа, — да, п. (Да. — да. 
где оа.=аТ/(46,) — значение параметра а в точке расщепления а. При 56,>2а 


(Ла=А%=Ао=0) петля (43) является разновидностью петли Люфт с Усилением без Усов 
(см. Рис. 206). При 6,<2а (лЛа=А=Аоз=0) петля (43) является разновидностью петли 
Люфт с Ослаблением без Усов (см. Рис. 20г). При 6,=2а (ла=А=Аоз=0) петля (43) вы- 
рождается в разновидность петли Люфт без Усов (см. Рис. 20е). 

На Рис. 196 изображена построенная по формулам (43) кусочно-линейная петля Лист 
при различных значениях фазового сдвига Аа/ (6,>2а, А=Ао%=0). Как видно из рисунка, с 
помощью фазового сдвига Ася возможно управлять усилением 7 Показанные на рисунке 


углы Ви уопределяются по формулам° 








Бад , _ Блапу _ Би, (©. + а.) + Ш, а. | 
Ь-а 6 -2а БШ (а +ла)+щ. а. Ш. Аа | в 
в Ь, (6, —2а) _ (6, -2а\апВ _ 6, [ых (а. + да )- ш. а. | | 
(6 -а) | Би. (м. + да) + 1. а. И. Аа | 


Используя формулы (45), уравнения (43) можно при необходимости выразить через углы 
Ви/или у Если задать фазовый сдвиг Ася равным -1.0 (7=41°), -0.73 (2=25°) и -0.5 (7=0°), то 
с помощью формул (43) можно построить такое же семейство петель с теми же парамет- 
рами а, 6», 6, (Бу<2а, А=Ас=0), что и на Рис. 19а, петли на котором построены по фор- 
мулам (40). Задавая ненулевые фазовые сдвиги Ао, До, Ао, с помощью уравнений (43) 
можно построить кусочно-линейные петли, показанные на Рис. 206, г-е, и, к. 

Перенося в формулах (43) параметр расщепления а в аргумент функции треугольных 
импульсов, можно получить систему уравнений, описывающую кусочно-линейную петлю 


гистерезиса Люфт с Усилением (даи=А%=А=0)° 
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Улучшенная параметрическая модель для аппроксимации петли гистерезиса 
х(<) =Ь 4. а, 


У(а)= (6, -Блапу п, Ел И, РН. А +6 фату, <. ее 
апр -—1апу 8 {апб-—1апу 8 


Параметры петли 6,, 6,, р, усвязаны между собой следующим образом: 26/6,=апр\|апу 
Указанная связь делает уравнения (46) существенно менее универсальными, чем давае- 
мые ниже уравнения (47) и (54). Тем не менее, с помощью уравнений (46) можно строить 
петли подобные петлям, изображённым на Рис. 20а-е. 
а. Универсальная формула кусочно-линеиной петли 

В ходе настоящего исследования получено общее выражение, описывающее кусочно- 
линейную петлю гистерезиса Люфт с Усилением"" (см. Рис. 20а)? 

х(<) =Б 4й. <, 
а. {апр 

{ап бр —{апу 


(47) 


у(<) = (6 -Ь лап у)ир. [с — | +6 лапу. а, 


в котором верхнее основание а трапецеидальных импульсов {гр. определяется по фор- 
муле 
_ Татр -Ь,) 
_ 26, (апр —1апу)’ и 
а нижнее основание О по формуле П= Т-4. Пошаговый вывод формул (47) дан в дополни- 
тельном материале. 
Кусочно-линейная петля (47) построена на комбинации треугольных и трапецеидаль- 
ных импульсов. При необходимости трапецеидальные импульсы {гр. в (47) можно пред- 


ставить в виде суммы двух треугольных импульсов® 


мы 
паи и 
Е ЖИ ИСИ га 
АЙ Е а И 
{ап б—1апу 2 {ап б—1апу 2 








Уравнения (47) являются универсальными, поскольку позволяют получить весь набор 
кусочно-линейных петель Люфт и Реле, показанных на Рис. 20. Например, полагая в (47), 
(48) 7=0°, получаем систему уравнений для построения петли гистерезиса типа Люфт"“ 
(см. Рис. 20д) 

х(<) =Б Аи. а, 
уе) =Ь тр. (а а) г 
в которой верхнее основание @трапецеидальных импульсов р. определяется по форму- 
ле 
_ Т(Ь 1апр -6,) ви 
26 {ап В 
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Полагая в (47), (48) 6=90°, получаем систему уравнений для построения петли гисте- 
резиса Неидеальное Реле с Усилением "° (см. Рис. 206) 
х(<) =Ь 4. <, 
у(<) = (6. -Ь лап удгеси. (и- а.) +6 лапу. а, (52) 
в которой верхнее 4 и нижнее О основания трапецеидальных импульсов одинаковы и 
равны 7/2. Последнее означает вырождение трапецеидальных импульсов в порождаю- 
щей функции у(а) в прямоугольные импульсы гес% со скважностью 2. 
Полагая в (52) 7=0°, получаем систему уравнений для построения петли гистерезиса 
Неидеальное Реле (см. Рис. 20и) 
х(<:) =Ь 4. ©, 
(а) =6 гесь, (и — а. ) (53) 
Для получения петель без усов (см. Рис. 206, г, е, х, и, к) один из параметров а, 6,, 6,, 
В в формулах (47), (50), (52), (53) должен выражаться через остальные параметры со- 


гласно соотношения 1апб=Ь/(5,-а). Для получения петель с ослаблением"” "8 


(СМ. 
Рис. 20, г, з, и) в формулах (47), (52) следует задавать отрицательный угол и 
На Рис. 21 даны примеры разложения кусочно-линейных петель Люфт с Усилением 


(6 х=а1 б/а, бч,у=а: б/а, рн = рб, я 


1 





й=х ал=а-аа, — Б2х=б-Ёх, 
Б2у=Бу-биу, б2=р, =, Люфт 
(б1х=а:буа, 
{апри= Та /[(Т-2а)61х|, и=0, 
ар=а-а1, 6>х=6х-блх, 
62,=бу-бчу, 

{апр>= ТБ2у//[( Т-24)62х|, }=0), 


Неидеальное Реле с Уси- 























лением (6:,=а.б/а, б=90°, 
И=Х  а2=а-а1,  6>х=бхЁчх, 
02,=бу-6- у б>=90°, = Я И 


Неидеальное Реле 
(6 х=а1 б/а, бн=90°, И=0, 
а›=а-а1, 6>х=бх-бчх, 


Ь2у=Бу-Вч, б2=90°, №=0) на - й | 

















0 0.5 1х(а) 





две петли того же типа, что 

Рис. 21. Разложение кусочно-линейной петли (зелёная) (а) Люфт 

с Усилением, (6) Люфт, (в) Неидеальное Реле с Усилением, (г) 

скольку каждую из петель Неидеальное Реле на две петли того же типа, что и разлагаемая 
петля. 


и разлагаемая петля.° По- 
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У(а) 
в 


У(а) 
1 = 


разложения в свою очередь 


можно разложить  анало- 


0.5 0.5 


гичным образом, то, как и 
гладкие петли, кусочно- 
линейные петли в общем 
случае можно представить 


в виде бесконечной суммы 














петель. 





На Рис. 22 даны приме- 





ры построения более слож- 
ных кусочно-линейных пе- 
тель: Люфт-Люфт-Люфт с 


0.5 





Ь1х = 0.30 




















| | 5-03 - Усилением"“ (р=2.4, 
В1= 90° 
ул = 0° 
|1 6: х=бмр, б1у=бу/р, ри=р, ПЕР 
р Ре а> = 0.28 — 
Бо» = 0.70 = =Ох- 
Ь2, - 085 ееЕЬ Веебк Их 
Во = 81° 
— Я РЬ о" - б2у=Бу-Влу, б2=р, р=р гдер- 
| - ы" 05 6 0.5 ' в общем случае произ- 
(в) (г) 


Рис. 22. Кусочно-линейные петли гистерезиса (а) Люфт-Люфт- вольное число), — Люфт- 
Люфт с Усилениаем, (6) Люфт-Люфт-Люфут, (в) Люфт-Реле-Люфт с 

Усилением, (г) Люфт-Реле-Люфт, полученные сложением двух Люфт-Люфт (р=2.4, 
петель типа (а), (6) Люфт, (в), (г) Реле и Люфт. Бир, у =БИр. В-В, м0, 
а›=а-аи, 6>х=6-6и» 6,=Бу, =В, р=0), Люфт-Реле-Люфт с Усилением” (6!=р, и=0, 
а›=а-а1, б>х=бх-бч», б2у=бу-бчу; а1, (> и } определяются численным решением системы 
уравнений°), Люфт-Реле-Люфт (2:=Д, и=0, а›=а-аи, 6>,=6-Ё», 6ьу=Бу-6у, 6=0; а! и Дь опре- 
деляются численным решением системы уравнений°). Петли типа Люфт-Люфт-Люфт 
представляют собой комбинацию из пары петель типа Люфт, петли типа Люфт-Реле- 
Люфт представляют собой комбинацию из петель типа Неидеальное Реле и типа Люфт.° 
Подробности образования/декомпозиции кусочно-линейных петель гистерезиса, пред- 
ставленных на Рис. 21 и Рис. 22, даны в дополнительном материале, там же можно найти 
разновидности этих петель с ослаблением и/или без усов. 

Обнаружение в ходе физических измерений петель, изображённых на Рис. 22, может 
служить указанием на то, что в рассматриваемой системе действует два разных нало- 
женных друг на друга гистерезисных процесса. 

р. Порождающая функция с пороговым элементом 
Помимо формул (47) для описания кусочно-линейной петли гистерезиса Люфт с Уси- 


лением (см. Рис. 20а) можно применить уравнения, в которых комбинация треугольных 


импульсов в порождающей функции у(а) пропущена через пороговый элемент НР 
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х(<) =6 и. а, 
у(<) = 2(6,—6 лап Ин] 


а1ап В, и-+Ь 1апу-6, 
(а-6,)4ап В +6 лапу 


1 54 
И. "| — ы +6 лапу. а. (54) 
Пороговый элемент Н; это функция реальной (неидеальной) ступеньки, задаваемая сле- 


дующим образом 
0, [< 0, 


н (0 = - 0, (55) 


| 
ЕЁ, 
где [& — “длительность фронта” реальной ступеньки. Параметр Ё определяется по пара- 
метрам петли а, 6,, 6, б, усогласно следующей формулы? 
_ 2(Ь, 1апу-Б,) (56) 
’ (а-6, ап +Ь Чапу’ 
и, наоборот, если зафиксировать значение длительности фронта ступеньки & то по этой 


формуле можно определить один из параметров, пере- 





численных выше. При #->0 (6-—90°) функция Н; вырожда- 
ется в идеальную ступеньку Н (функцию Хевисайда). В 
этом случае образуются петли типа Неидеальное Реле 
(см. Рис. 206-к). о 
Как правило, уравнения (54) позволяют получить весь 


—0.5 


набор типов кусочно-линейных петель Люфт и Реле, по- 





казанных на Рис. 20. Однако, при некоторых значениях 


— 








параметров а, 6, 6, уна петле могут возникать дополни- - 5.5 о . 05 1 х(а) 
Е а РУ у - Рис. 23. Кусочно-линейная петля 





тельные изломы (см. Рис. 23, сравни с Рис. 20а), превра- 
щающие петлю Люфт в петлю Люфт-Люфт"°. 
При необходимости формулы (40)-(47), (50), (52)-(54) 


можно выразить через трапецеидальные импульсы, ис- 


гистерезиса Люфт-Люфт с Уси- 
лением. Дополнительные изло- 
мы возникают, когда расщепле- 
ние а и наклон Дб становятся 
меньше определенной величи- 
НЫ. 


пользуя следующее представление треугольных импульсов в виде суммы двух трапецеи- 


дальных (0=За, Т=а+0) 


#1, (@) = 5 Ир. (а Е я) о Бр. (а г 


В 


#1. (@) = Ир. (а+ а) + Ир. (а-— р 


в 


3. Сдвинутые петли 


В общем случае сдвинутые кусочно-линейные петли (см. Рис. 24) получаются сложе- 


нием двух кусочно-линейных петель (47). Параметры расщепления а!, а> в таких петлях 


задают горизонтальный размер и/ сдвинутого участка, параметр насыщения 6, опреде- 
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У(а) 
1 







0.5 




























Ь1х = 0.5 
1 = 0.54 = 
Вл = 90° 
71 = 17° 


й = 0.2 


а = 0.25 4. 
Бох = 0.5 
62, = 0.46 
Вг = 90° 
у2 = 17° = 

















4 0.5 0 0.5 1х(@) 


Рис. 24. Сдвинутые кусочно-линейные петли гистерезиса (а) 
Люфт с Усилением, (6) Люфт, (в) Неидеальное Реле с Усилени- 
ем, (г) Неидеальное Реле, полученные сложением двух петель 
типа (а), (6) Люфт, (в), (г) Реле. 


р ии | {апу 
а=а, — а 1, 
{апр 2\1апр 


а расщепление а! — по формуле 








а. = 
| 21ап р, 


[а _ о (ал В, —1апу)+ Б, 1апу-26,, 


и гистерезиса 


ляет вертикальное положе- 
ние ЛР этого участка, усиле- 
ние уустанавливает наклон 
участка, а угол 6! (62=р!) 
определяет такой наклон 
складываемых петель, что- 
бы общий наклон В’ сум- 
марной петли оказался ра- 
вен заданному. 

Для петли Сдвинутый 


19, 20 ( 


Люфт с Усилением СМ. 


Рис. 24а) расщепление а 


находится по формуле 


(58) 


(59) 


Остальные параметры а», В, д! определяются численным решением следующей систе- 


мы уравнений (61„=6,/2, и=у 6>х=Б-Ьч», 6у=Бу-Ьну, =, в=9 


=тап 2, 


где х(а), у(а) — порождающие функции сдвинутой петли; @, , 


(60) 


@ — значения параметра & 


а 


в точках а\, а» сдвинутой петли, соответственно; с. (а,), ©. (В.), а. (В,) — зависимости 
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параметра а в точках еф, 6 
сдвинутой петли от пере- 
менных а2 и ру, соответст- 
венно; У. „6, — зависи- 
мость порождающей функ- 
ции от ББ, в точке ар 


(а=а, ).” 

















На Рис. 246, в, г приве- 
дены примеры сдвинутых 
петель Люфт,“ Неидеаль- 


ное Реле с Усилением„”' 


2 

а. = -0.10 
Неидеальное Реле,” яв- т, 
у“ 


В1= 90° 
ляющихся частными слу- 17° 


чаями петли Сдвинутый > и 
77$ О 


Б.,, = 0.30 
Люфт с Усилением. Другие о, 
ур = 17° - 
разновидности сдвинутых 


—1 -0.5 0 0.5 1Х( <) -1 -0.5 0 0.5 1Х( 




















петель с ослаблением ри’ ок, Трёхуровневые кусочно-линейные петли гистерезиса (а) 


и/или без усов даны в до- Люфт с Усилением, (6) Люфт, (в) Неидеальное Реле с Усилени- 
ем, (г) Неидеальное Реле, полученные сложением двух петель 
полнительном материале. типа (а), (6) Люфут, (в), (г) Реле. 


Обнаружение в ходе физи- 
ческих измерений сдвинутой петли указывает на то, что в рассматриваемой системе, по- 
видимому, одновременно действует два отдельных наложенных друг на друга гистере- 
зисных процесса (например, имеется плёнка из двух слоёв, магнитные свойст- 
ва/состояния которых различаются"! =“). 
а. Трехуровневые петли 

Частным случаем сдвинутых петель являются трехуровневые петли (см. Рис. 25). 
Трехуровневые петли возникают при задании отрицательного расщепления а и выполне- 
нии следующего равенства° 


2а{ап рб 
ай (61) 
{апу — {ап в 
Описание треёхуровневых петель с ослаблением и/или без усов имеется в дополнитель- 
ном материале. Обнаружение в ходе физических измерений трехуровневой петли может 
являться признаком того, что в рассматриваемой системе одновременно действует два 


различных наложенных друг на друга гистерезисных процесса. 
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С. Двойные петли гистерезиса 
1. Сцепление петель в точке насыщения 

Для более адекватного описания двойных петель гистерезиса способ, предложенный 
в работе 1, был усовершенствован. В частности, уравнения двойной гладкой петли само- 
непересекающейся в точке начала координат (петля в форме 0) согласно усовершенст- 


вованному способу выглядят следующим образом (а=0...2п): 


х(а) = $ 2а- зол —@)- ло) +6, здп(л — а) = (2тефа- \2а- да. — = + ь 
(62) 
у(@) = га воли -а)- ло) + 6, здп(т — <) = (тема - {Уго да, — ы + ь 


где з9апа=а/|Ч -— сигнум-функция; гео=Н(а-Н(а-л) — прямоугольный импульс шириной д. 
Двойная петля (62) образуется путём сцепления двух петель в точке насыщения б56, где 
порождающая функция у(@) достигает максимума 6,. На Рис. 26а показан пример двойной 
гладкой самонепересекающейся петли, построенной по формулам (62). 

Согласно (62) движение по двойной петле начинается в точке (0, 0) и происходит про- 
тив часовой стрелки сначала по верхней петле затем по нижней, движение заканчивается 
в точке (0, 0). Отличие двойной петли (62) от использовавшейся ранее заключается в том, 
что точка начала второй петли совпадает с точкой окончания первой, а точка окончания 
второй с точкой начала первой. Гаким образом, двойная петля (62) при изменении пара- 
метра хот 0 до 2л прочерчивается непрерывно (“не отрывая карандаш от бумаги”). 

Петля, при движении по которой происходит самопересечение в точке начала коорди- 


нат (петля в форме 8), строится по формулам (а=0...2т): 


х(а) = 4 [22 й ы зап(л - &)- ль + 6, зап(л — а) 
= гес1 “2 —Аа; — = + ь, + (1-—гес “| Аа; — 2а = ь 
у(а) = ; [22 ы здп(л — @)- ла + 6, здп(л — а) 


= гес! 4 У22- Аа. — = + ь, + (1 гес! «= — да. — 2а — ь 


Внешне петли (62) и (63) не отличаются друг от друга. 


(63) 


Двойные самонепересекающиеся и самопересекающиеся кусочно-линейные петли 
строятся по формулам (62) и (63), соответственно, в которых л заменён на Г/2 (0=0...Т). 
На Рис. 266 показан пример двойной кусочно-линейной самонепересекающейся пет- 


16, 19, 24 
) 


ли построенной по формулам (62) при Ааз=0. Самопересекающаяся петля внешне 


выглядит также. 
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(а) (б) Рис. 27. Двойная гладкая петля 
Рис. 26. Двойная (а) гладкая, (6) кусочно-линейная петля гисте- гистерезиса, образованная сце- 
резиса, образованная сцеплением в точке насыщения 6 двух пе- плением двух петель Классиче- 
тель типа (а) Классическая, (6) Люфт с Усилением. В точке сцеп- ская в точке х„.„. В точке сцеп- 
ления петлю можно сделать как самонепересекающеися, так И ления петлю можно сделать как 
самопересекающейся. самонепересекающейся, так и 
самопересекающейся. 








2.Сцепление петель в точке Хтах 

Вместо сцепления двух петель в точке насыщения 6, петли можно сцеплять в точке, 
где порождающая функция х(а) достигает максимума хт„». Чтобы определить значение 
параметра ош.» при котором х(@атах)=Хтах следует решить (численно) уравнение 


ах(<)/ а = 0 а. ИЛИ в развернутом виде 


тазт(а „+ Аа )соз” („+ да)- пб, с0$(и „+ Аа) т" "(+ А) = 0. (64) 


Двойные самонепересекающиеся петли, сцепленные в точке хи„х, строятся по формулам 


х(а) = {2а- 5 (69п(т - а) +1) + ра +х(а„)здп(л -а) 


= (1- 2'ема)(ха+а„)- ха} 


у(а) -- У 2а- 5 (69п(т ыы &)+ 1) + ны, т У(а„„х }9п(л — а) 


= (1- 2ема)(у(2а-+а_)-у(а„,)) 


тах 


Двойные самопересекающиеся петли, сцепленные в точке хтах, строятся по формулам 


х(а) = 4 [22 _ ы зап(л-а)-+а- : + ха )зап(л - а) 
))— геста(х(2а+ а „)- ха), 


тах 


= (1-геча) (жа, -2а)- х(а 


тах тах 


. (66) 
у(@) — ; [22 = зап(л —_ а) т тах — ы Е У(а„„х}здп(л — а) 
-- (1 — гес а)(у (ох — 24) _У (ах )) — гес (у (20 я так) У (ах )) 
На Рис. 27 показана двойная самонепересекающаяся петля, построенная по формулам 


(65). Двойная самопересекающаяся петля, построенная по формулам (66), внешне вы- 


глядит также. 
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Рис. 28. Двойная самонепересе- (а) 





0.5 1 х( а) 


кающаяся петля гистерезиса РИС. 29. Двойные самонепересекающиеся петли гистерезиса ти- 
типа Пропеллер, образованная па Пропеллер, образованные передавливанием петли в точке 
путём замены горизонтального Начала координат посредством (а) возведения в степень, (6) ну- 


расщепления на вертикальное. Левого расщепления а и фазового сдвига Л4. 


3. Замена горизонтального расщепления на 
вертикальное 
Простейший способ получения двойных самонепересекающихся петель типа 
Пропеллер"? основан на представлении (2), в котором вместо расщепления а (а=0) по х 
выполняется расщепление а, по у. Уравнения расщепляющей кривой в этом случае 


имеют следующий вид: х2(а=асоз”а у2(а=аусоз”а Из уравнения У а, , )=Ву2, 


записанного для полувысоты точки насыщения Б, (вместо полувысоты можно 
использовать любое нужное значение), находим зависимости (т=1) 


Ь, Е + \/4а, +36, (67) 


а, (а,)= агссо$ 2 
у у 


После чего, решая численно уравнение 
ха, (а, )- Ма. (а, = 2а., (68) 

определяем такое значение вертикального расщепления а, при котором величина гори- 
зонтального расщепления аз двойной петли на полувысоте точки насыщения 5, будет 
равна заданному значению. Полученная описанным способом двойная петля Пропеллер 
показана на Рис. 28. Рабочие выражения для случая т-3, 5,... даны в дополнительном 
материале. 
4.Передавливание петли в точке начала координат 

Двойные самонепересекающиеся петли типа Пропеллер можно образовывать “пере- 
давливанием” петли в точке начала координат (см. Рис. 29а). Передавливание происхо- 
дит при возведении порождающих функций х(а), У(а в нечётные степени К и / согласно 


выражениям (27). Из уравнения у|а.)=6,/2, записанного для полувысоты точки насы- 
а у 


щения 5, (вместо полувысоты можно использовать любое нужное значение), находим па- 


раметр а 
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| 1 
о, = асзт———. 
а | 2ВЬ!" (69) 
Решая уравнение 

Х(а,}- х(л - а,)= 2а,, (70) 
можно определить такое значение расщепления а, при котором величина расщепления аз 
двойной петли на полувысоте точки насыщения 6, будет равна заданному значению. 


Аналитическое решение уравнения (70) для К=З выглядит следующим образом 


2 
у 2АВ, + \/4А7 6} тб” ©, + 22} —- 2 АБ; т?” ©, (71) 
с0$” 314 А?(а, + \/4 АБ? зп?” и, + а: | 


На передавливании петли основан ещё один способ получения двойных самонепере- 
секающихся петель. В этом способе расщепление а в формуле (13) задаётся равным ну- 
лю, а фазовый сдвиг ДАаф (или Дов), устанавливается отличным от нуля. В таких петлях 
сдвиг Ао (или А0з) выполняет функцию расщепления а. На Рис. 296 показан пример 


двойной самонепересекающейся петли типа Пропеллер,” 28 


полученной передавливани- 
ем. Из уравнения У(а)=6,„2, записанного для полувысоты точки насыщения 6, (вместо 
полувысоты можно использовать любое нужное значение), находим параметр а=л/б. 
Решая численно уравнение 

х(а.)- х(т- а.) = 2а,, (72) 
можно определить такое значение фазового сдвига Лор (или Аоз), при котором величина 


расщепления аз двойной петли на полувысоте точки насыщения 6, будет равна заданно- 





му значению. 


У [72 
1 


5. Двойная кусочно- 
линейная петля как м 
частный случаи кусочно- 
линейной сдвинутой о 
петли 

0.5 


Задавая нулевое рас- 


щепление а=0, сдвинутая 














петля (см. раздел И.В.З) 





(а) 

вырождается в двойную Рис. 30. Двойные —кусочно-линейные —самонепересекающиеся 
петли гистерезиса Люфт с Усилением (а) с перемычкой, (6) без 

самонепересекающуюся перемычки. Петли образованы в результате сложения двух оди- 

петлю Люфт с Усилением?” Нарных петель Люфт с Усилением. 

(см. Рис. 30а). Задавая помимо нулевого расщепления высоту сдвинутого участка соглас- 


но формулы Р=и{апУ2, можно получить двойную самонепересекающуюся петлю Люфт с 
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Усилением без перемычки“ ^^ 


(см. Рис. 3060). Двойные самонепересекающиеся петли 
Люфут, Реле с Усилением и Реле, имеющие и не имеющие перемычку, а также разновид- 
ности этих петель с ослаблением и/или без усов даны в дополнительном материале. Из 
способа образования двойной петли следует, что её обнаружение в ходе физических из- 
мерений может являться признаком того, что в рассматриваемой системе одновременно 


действует два разных наложенных друг на друга гистерезисных процесса. 


О. Троийные петли гистерезиса 
1. Сцепление петель в точках насыщения 
Тройную петлю можно собрать из трех петель — центральной и двух внешних. Трой- 


ные гладкие петли, сцепленные в точках насыщения 56, строятся по формулам (0=0...2х) 


х(е) = (ема-тесКа-л))х, [за _ = 


+ (‘ест а-+тесКа-л)- Ч х |+ За-— = (6 +6, )зап(л -— “) 
(73) 


у(а) = (еска+тесК(а-л)у. [за- Е | 


+ (геста+тесКа-—л)-— Чу. и За-— ы _ (6, +6,, )здп(т — “) 


где х!(а, у((а) - уравнения центральной петли; х2(а), у›(а) — уравнения внешних петель; 
бчх, бл, — координаты точки насыщения центральной петли; 62», 62, — координаты точек на- 
сыщения внешних петель; гео=Н(а)-Н(о-л/З) — прямоугольный импульс шириной л/З. 

При сборке петель (73) обычно выполняется условие и=»ь, где и, 75 — углы наклона 
касательных к нерасщепленной (а.=0) центральной петле и нерасщепленным (а>=0) 
внешним петлям, соответственно, в точке насыщения 61. Углы наклона и, 6 касательных 
определяются по формулам° 


Ь, 
у. = агап—^—, 
Ш 1х 


(74) 


о) 
у, = амап—^”_. 
П> 2х 


Если в уравнениях (73) аргумент За функций х2(а), у›(а) используется со знаком плюс, 
то получается самонепересекающаяся петля, если со знаком минус — самопересекаю- 
щаяся. Внешне самонепересекающиеся и самопересекающиеся петли друг от друга не 
отличаются. Типы центральной и внешних петель могут быть различными. На Рис. З1а 
показана тройная гладкая самонепересекающаяся петля, построенная по уравнениям 
(73). Петля состоит из центральной Классической петли (1) (а!=а) и пары внешних Клас- 


сических петель (1) (6>,=(Б»-6чх)/2, б>у=П»б>4апи, ь=и). В дополнительном материале да- 
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ются формулы тройных пе- 
тель, сцепленных в точках 
Хтах. 

Тройные самонепересе- 
кающиеся и самопересе- 
кающиеся кусочно- 


линейные петли’? получа- 

















ются из формул (73) заме- 


а) — | | | (6) — | | | ной л на 7/2 (0=0...Т). На 
Рис. 31. Тройная (а) гладкая, (0) кусочно-линейная петля гисте- 


резиса, образованная сцеплением в точках насыщения 6 трёх Рис. 316 показана тройная 
петель типа (а) Классическая, (6) Люфт с Усилением. В точках 
сцепления петлю можно сделать как самонепересекающейся, так 
и самопересекающейся. пересекающаяся петля, по- 





кусочно-линейная самоне- 








строенная по уравнениям 
(73). Петля состоит из цен- 
тральной петли Люфт с 


Усилением (47) (а!=а) и па- 





ры внешних петель Люфт с 
Усилением (47) 
(=), Бу=(Ьу-Вну)/2, 
№=И). 

2.Петли с усами 























Рис. 32. Симуляция одинарной Классической (а) гладкой, (0) гиб- 
ридной петли с длинными усами с помощью тройной самонепе- произвольно большой 
ресекающейся петли гистерезиса. Усы — внешняя пара нерасще- 
плённых петель типа Лист (а) гладких, (6) кусочно-линейных. 


длины 

Тройные петли гистере- 
зиса (73) удобно применять, в тех случаях, когда требуется получить одинарные глад- 
кие. 3 или одинарные гибридные петли гистерезиса с длинными усами. Напомним, что 
длинные усы в модели (1) можно получить, увеличивая степень т,° однако одновременно 
с этим происходит значительное изменение кривизны петли. На Рис. З2а показана симу- 
ляция одинарной гладкой петли наклонная Классическая с длинными усами с помощью 
тройной самонепересекающейся гладкой петли (73). Гройная петля состоит из централь- 
ной наклонной Классической петли (30) (а!=а, 9=0) и усов, образованных из пары внеш- 
них нерасщеплеённых петель типа Лист (1) (а2=0, 6>х=(6х-6чх)/2, Бу=(Бу-б1у)/2, Тр=З, Пг=1), 
ориентированных под углом 7=агап(6>/6>х). 

Необходимое искривление к: петли определяется по формуле 
Ь,[апб,(1-п, 6, -В,)-Б,+6,|+16,(6,-Ь,) 


к, = агЧап 
Ь,, (п, - Тапа (6, -6,)-6,+Ь, 


(75) 
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Формула (75) найдена из условия и=?, где и, 75 — углы наклона касательных к нерасщеп- 
лёенной (а.=0) центральной петле и нерасщепленным (а2=0) внешним петлям, соответст- 
венно, в точке насыщения 61. Из формулы (75) для некоторого произвольного значения х: 
можно определить соответствующее значение 6 или для произвольных значений 6,, ж, 
9; найти соответствующее значение 5,. 

На Рис. 326 показан пример построения из трёх петель наклонной гибридной Класси- 
ческой петли с усами.” "? Петля состоит из центральной наклонной гибридной Классиче- 
ской петли без усов (38) (а!=а, 61,=бу-262,, 9=0, и=7, к=к) и усов, образованных из пары 
внешних нерасщеплеённых кусочно-линейных петель Люфт с Усилением без Усов (40) 
(а›=0, 6>х=(бх-бчх)/2, б>у=б2апь, Тр=Пр2=1, ь=и).5 В связи с переносом в (31) точки насы- 
щения 6 из позиции ад= 1/4 в позицию а= 1/8 (Аа= 1/8) в формулах (73) помимо замены л на 
Т/2 из аргументов функций х! и у! следует вычесть величину 1/8. Кусочно-линейные петли 
с усами любой длины строятся непосредственно по уравнениям (47). 

3. Добавление вертикального расщепления 

Простейший способ получения тройных самопересекающихся петель типа Классиче- 
ская основан на представлении (2), в котором помимо расщепления а, по х выполняется 
дополнительное расщепление а, по у. Уравнения расщепляющей кривой в этом случае 
имеют следующий вид: х2(а)=ахсоз”а, у›(@=аусоз”а. Уравнения тройной петли с учётом 
дополнительного расщепления записываются как 


х(<) = а, соз" х+Б зп” @, 


‹. Ре (76) 
У(а) = а, со$” аб эта, 
где параметры а», а,, 6,, О, определяются по формулам (т=1):8 
_ асо$” а, -6 9” <. 
о эм” а, +с0$"" а, 
а =-6. $зпа 
у у а’ 
(77) 


азт а. +6,.С0$ 4. 
эт” и, + 0$" а. 


с > 


х 
Ь, =6,С0$ @., 

где а. - значение параметра ав точке расщепления а. Отрицательное с. задает величи- 

ну перехлеста тройной петли. Рабочие выражения для случая т=3, 5, ... даны в дополни- 

тельном материале. Как следует из формул (77), при а.=0 тройная петля вырождается в 

обычную петлю Классическая. Задавая а.>0, можно изменять форму и кривизну одинар- 

ной петли.? 


25, 29 


Полученная описанным способом тройная петля Классического типа показана на 


Рис. 33. Стоит отметить, что, полагая расщепление а=0, тройная самопересекающаяся 
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Рис. 33. Тройная гладкая само- 
пересекающаяся петля гистере- 
зиса Классического типа, обра- 
зованная в результате дополни- 
тельного вертикального расще- 
пления. 


ную самонепересекающуюся 


4. Сжатие с перехлестом 





—0.5 


— | 
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Рис. 34. Тройная гладкая само- 
пересекающаяся петля гистере- 
зиса Классического типа, обра- 
зованная в результате “сжатия” 


фазовым сдвигом Д(ь. 


петля вырождается в двой- 


Р. В. Лапшин 








а1 = -0.12 
Бух = 0.5 
Ь1у = 0.67 -— 
В = 83° 

71 = 17° 


а = 0.25 — 
Бох = 0.5 
62, =0.33 
До = 83° 
2 = т" — 








| | | | 
0 0.5 





= —0.5 1Х( а) 
Рис. 35. Гройная кусочно- 
линейная самопересекающаяся 
петля гистерезиса типа Люфт с 
Усилением, образованная в ре- 
зультате сложения двух одинар- 
ных петель Люфт с Усилением. 


петлю типа Пропеллер (см. также раздел И.С.3).? 


Тройные самопересекающиеся петли в улучшенной модели (13) образуются путем 


задания отрицательного фазового сдвига А (или положительного До), настолько 
сильно “сжимающего” петлю, что возникает перехлест. Пример тройной петли Классиче- 


^>, = полученной данным способом показан на Рис. 34. Стоит отметить, что, 


ского типа, 
полагая расщепление а=0, тройная самопересекающаяся петля вырождается в двойную 
самонепересекающуюся петлю типа Пропеллер (см. раздел 1.С.4). 
5. Тройная кусочно-линейная петля как частный случай кусочно-линеиной 
сдвинутой петли 

Из сдвинутой петли (см. раздел ||.В.3) можно получить тройную самопересекающуюся 
петлю (см. Рис. 35). Для этого расщепление а должно быть отрицательным, а горизон- 
тальный размер и/ сдвинутого участка петли должен превосходить значение их, опреде- 
ляемое формулой (61). Тройные самопересекающиеся петли Люфт, Реле с Усилением, 
Реле, а также разновидности этих петель с ослаблением и/или без усов даны в дополни- 
тельном материале. Из способа образования тройной петли следует, что её обнаружение 
в ходе физических измерений может являться признаком того, что в рассматриваемой 
системе одновременно действует два разных наложенных друг на друга гистерезисных 


процесса. 


Е.Площадь петли гистерезиса 


Площадь петли гистерезиса характеризует потери энергии в пьезоэлектри- 


ке/ферромагнетике при приложении к нему переменного электрического/магнитного поля. 


Чтобы найти площадь петли гистерезиса воспользуемся известной общей формулой 
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И 


Согласно (78), масштабирование х(а)= Ах(а), у(а) = Ву(а) петли гистерезиса приводит к 
изменению её площади пропорционально произведению масштабных коэффициентов Аи 
В, т.е. $ = АВ$. На Рис. Зба показана серия из трёх петель гистерезиса, площади кото- 
рых при масштабировании удваиваются. 

1.Гладкие петли 


Подставляя в (78) выражения х(а), у(а) улучшенной модели (13) вместе с их произ- 








водными и интегрируя, получаем” °° 
^ т-1 Г п-1 
$=| ЧС 2 соз(ла, — лм) +—*С 2 эт(ла — д) ль 
рт-1 1 3 рп п 2 3 у 


| (79) 
В О У 
\Мап” ф,/+1 


где п — нечётное число; А1, В+ — коэффициенты Фурье для первой гармоники (см. форму- 
лы (25)). Из формулы (79) следует, что площадь петли (13) определяется амплитудой Ат 
и фазой $ только первой гармоники; остальные гармоники разложения (7) порождающей 
функции х(а) никакого влияния на площадь $ не оказывают. 

Наклон петли с помощью фазового сдвига Ао (т>1, А%=Аоз=0) вызывает рост её 
площади $ (см. Рис. 3). Площадь $ петли гистерезиса возрастает при увеличении фазо- 
вого сдвига Лор (см. Рис. 4). При увеличении фазового сдвига Ао площадь $ уменьшает- 
ся (см. Рис. 5). 


Из формулы (79) следует, что при А=Ао площадь петли 
т 
ся С лав, | (80) 

2” соз" "(да — Да, ) 
Таким образом, при указанном условии петли любого типа (пл не входит в формулу) с лю- 
быми значениями насыщения 56, будут иметь одну и ту же площадь $ при одних и тех же 
параметрах а, бу и т (см. Рис. 366). 

Из формулы (80) следует, что в случае А=Лаз=0 площадь петли с отрицательным 
наклоном (Аси>0) в точке расщепления равна площади петли с положительным наклоном 


(Ла.<0) в точке расщепления, т.е. $. — . Из формулы (80) также следует, что 


а 


З Таким образом, для любого типа петли в случае д%=Аоз=0 (т>1) 
1 


-лаз >0 ы а | 
потери энергии в наклонённых петлях (Аои=0, см. Рис. 3) больше чем в “прямых” (Ади=0). 


Из формулы (80) также следует, что при Т=1 (А=А оз) площадь петли 
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$ = лаб... (81) 


Таким образом, при указанных условиях петли любого типа (п не входит в формулу) с 
любыми значениями насыщения 6, и с любыми значениями фазовых сдвигов будут иметь 


одну и ту же площадь при 





_ у(а) 
1 


3 
_— 
и. 


одних и тех же параметрах 


аи 6,. В качестве примера 


0.5 
0.5 


на Рис. З6в представлены 
петли одинаковой площади ° 
с разными фазовыми сдви- о+.- 


—0.5 


гами Ао. Формула для оп- 








ределения площадей пе- 




















тель Лист (п=1), у которых (а). 








т=1, а Ла=Ао, выглядит 
также как и (81). Площади 
этих петель также не зави- т 
сят ни от величин насыще- 
ния б», ни от величин фазо- 


вых сдвигов. Площадь (81) 





численно равна площади 














| | | | | 
1 х( а) —1 —0.5 0 0.5 1 х( <) 


эллипса с полуосями аи 6,. 











(г) 

Рис. 36. Площадь петли гистерезиса (а) пропорциональна произ- 

щади петли, построенной В@едению масштабных коэффициентов А и В по осям хи у, соот- 
ветственно, (6) одна и та же при любых насыщениях 6, (ограни- 

по модели (1), в формулу чивающее условие: А=Л ов), (в) одна и та же при любых фазо- 

вых сдвигах Аа (ограничивающие условия: т=1, А%=Аов), (г) 

одна и та же при любых насыщениях 6,. 


Для вычисления пПлЛоО- 


(79) следует подставить 
нулевые фазовые сдвиги 


Ада=А0=Аоз=0. В результате получаем 





$= С„? лаб, = ——=^— =лАф,. (82) 


В работе 1 в формуле (27) и в сопроводительном тексте к этой формуле допущена 
неточность. Формула (27) даёт верный результат для любого п — чётного или нечётного. 
Поскольку формула (27) не даёт решение при т=1, то вместо формулы (27) следует ис- 
пользовать формулу (82), приведённую в настоящей работе. По той же причине вместо 
формулы (28) работы 1 следует использовать соответствующую формулу, приведенную в 
дополнительном материале. 

Поскольку п не входит в формулу (82), то у всех трех типов петель Лист, Месяц и 


Классическая модели (1) площадь петли $ — одинаковая при условии, что остальные па- 
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раметры а, Б, и т этих петель одни и те же (см. Рис. 1). Поскольку в формулу (82) пара- 
метр 6, также не входит, то все петли (1) с одинаковыми параметрами а, 6, и т, но раз- 
ными 6, имеют одну и ту же площадь $ (см. Рис. З6г). 

По формуле (82) также вычисляется площадь гладких петель (29), наклоненных в точ- 
ке расщепления на угол 0 посредством перекоса системы координат. Поскольку угол пе- 
рекоса 9в формулу (82) не входит, то петли, наклонённые перекосом на любой угол д, у 
которых параметры а, Б,и т одни и те же, имеют одну и ту же площадь $ (см. Рис. 8). 

Площадь Классической петли (30), наклоненной и искривлённой перекосом, вычисля- 
ется по общей формуле? 

ты 
7  2К+1 


ла 
э = б {ап 1-— (т+1) +6 +, 
а тик] (83) 





где п — нечётное число. Согласно (83) площадь не зависит от угла перекоса 0, поэтому, 
петли, наклонённые перекосом на любой угол 6, у которых все остальные параметры од- 
ни и те же, имеют одинаковые площади (см. Рис. 9). При хк=0 (см. Рис. 8) формула (83) 
трансформируется в (82). В частном случае, например, при т=п=3З (см. Рис. 9) формула 


(83) выглядит следующим образом 
$ = Е ла(Ь, {апк +26, ) (84) 


Интересно отметить, что при Ла=А%=Ааз=А 00, где Ах - произвольное действитель- 
ное число, формула (79) также приобретает вид (82). В этом случае петли, построенные 
согласно моделей (1) и (13), имеют одинаковые площади. Дело в том, что при указанных 
условиях эти петли совпадают по форме друг с другом (см. раздел П.А.2.а). 

Площадь Классической петли (14), наклонённой поворотом, вычисляется по общей 
формуле? 

Е 
т 2 2К+1 


7: 
а п 96, с0$ 9 -— Ь, $п9]11- (т+1) +6 +, 
2° | вот+п-2к| И 





где п — нечётное число. Согласно (85) площади классических петель различны для углов 
поворота 6, имеющих противоположные знаки. При 9=0 формула (85) трансформируется 
в (82). В частном случае, например, при т=п=З формула (85) выглядит следующим обра- 


ОМ 
© = Бла б(Ь, 059 —6, 5п 9)+ 25, | (86) 
2.Кусочно-линейные и гибридные петли 


Площади простейших кусочно-линейных петель, имеющих форму параллелограмма 
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(см. Рис. 15, Рис. 19, Рис. 20, Рис. 21), проще всего определить, используя формулу пло- 
щади параллелограмма (т=п=1)° 


(а-Ь, 4апу+Ь, И ВЕЕЦУЕЩИ 


2—4 {ап —1апу р (8-97) | (87) 


При отсутствии усиления (7=0) формула (87) существенно упрощается 
5 = 4аб.. (88) 

Поскольку площадь (88) соответствует площади параллелограмма со стороной 2а и вы- 
сотой 26,, то по этой формуле можно определить площади кусочно-линейных петель типа 
Люфт и Неидеальное Реле (см. Рис. 20Дд, е, и, к). Поскольку гибридные петли (37) с нуле- 
выми фазовыми сдвигами (см. Рис. 11, Рис. 17) обладают трансляционной симметрией, 
то их площадь равна площади параллелограмма со стороной 2а и высотой 26, и поэтому 
также может быть найдена при помощи формулы (88). 

Площадь гибридной Классической петли с усилением/ослаблением (70) (38) (см. 
Рис. 18) вычисляется по следующей общей формуле? 


да ВИТ И ] + (6, - 8 = ‘апу—пап -). (89) 
) п+К 


$= ка -а/апк- 6, 1апу+ О 


где п — нечётное число, К — чётное. При отсутствии усиления (7=0) формула (89) обретает 
простой вид (88). Как следует из (88), площади гибридных петель с нулевым усилением 
не зависят от 6, п, дик. 

Поскольку значения параметра а во всех угловых точках кусочно-линейных петель из- 
вестны, то известны и декартовы координаты этих точек. Следовательно, площадь пе- 
тель-многоугольников, таких как Люфт-Реле-Люфт, Люфт-Люфт, Люфт-Реле (см. Рис. 12- 
14, Рис. 16, Рис. 23); Люфт-Люфт-Люфт, Люфт-Реле-Люфт (см. Рис. 22); Сдвинутый 
Люфт, Сдвинутое Неидеальное Реле (см. Рис. 24) и других им подобных может быть вы- 
числена по формуле площади многоугольника, заданного координатами своих вершин. 

В ряде частных случаев петли-многоугольники вырождаются в такие многоугольники, 
которые можно представить в виде совокупности трапеций и/или параллелограм- 
мов/треугольников. Например, петлю Люфт-Реле-Люфт с Усилением на Рис. 22в можно 
представить в виде двух одинаковых трапеций и прямоугольника, а петлю Сдвинутый 
Люфт с Усилением на Рис. 24а - в виде трёх параллелограммов, два из которых одина- 
ковые. Таким образом, для вычисления площадей подобных петель можно использовать 


соответствующие формулы площадей трапеций, параллелограммов, треугольников. 


Ш. ПРИМЕНЕНИЕ УЛУЧШЕННОЙ МОДЕЛИ 
С помощью разработанной модели можно строить гладкие, кусочно-линейные, гиб- 


ридные, частные, зеркально-отражённые, обратные, реверсивные, двойные и тройные 
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Таблица. Усредненная относительная 


аппроксимации (6) (%) исходной и улучшенной пара- 
метрических моделей петли гистерезиса. 


Петля 
Лист, 
Рис. 37а 


Петля Клас- 
сическая, 
Рис. 376 


погрешность 


Петля Клас- 
сическая, 
Рис. 37в 


петли. Вычисление производных, 


нахождение Гармонически линеари- 


функции 
гистерезисного звена и обратной 


зованной передаточной 


функции Производятся в улучшенной 


модели аналогично тому, как это де- 


лалось в работе 1. 





С ЕСИ ОКА 


Ас, А, Аз и ряда других преобразований, изменяющих наклон и кривизну петель, по- 


Использование фазовых сдвигов 


зволяет в несколько раз уменьшить погрешность аппроксимации петель гистерезиса. 
Анализ погрешности, проведённый по методике, описанной в работе 1, показал, что ус- 


реднённая относительная погрешность аппроксимации 
100% х у (х) 
(5) = ль 2" х 


216, 
— данные модели; у-(х) 








(90) 


(где уз(х) — данные эксперимента; п — количество точек на восхо- 
дящем или нисходящем участке экспериментальной петли) улучшенной модели не пре- 
вышает 1% (см. Таблицу).? Для корректности сравнения при определении погрешностей 
использовались те же экспериментальные петли гистерезиса, что и в работе 1. На Рис. 37 
показаны экспериментальные петли из работы 1, на которые наложены аппроксимирую- 
щие петли, построенные по существующей (1) и улучшенной (13) моделям. Легко видеть, 
что точность аппроксимации улучшенной модели заметно выше точности аппроксимации 
существующей модели. Для минимизации погрешности аппроксимирующую петлю реко- 
мендуется проводить, используя метод наименьших квадратов.?" 

Также как и исходная модель, улучшенная модель применима для исправления иска- 


жений, вносимых гистерезисом пьезоманипуляторов (см. Рис. 37а) сканирующего зондо- 
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Рис. 37. Аппроксимация реальных гладких петель гистерезиса типа (а) Лист, (6), (в) Классическая. 
Красная негладкая петля — экспериментальная, синяя пунктирная петля — существующая модель, 


зеленая гладкая петля — предложенная улучшенная модель. Погрешность аппроксимации (5) 
улучшенной модели не превышает 1%. 
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вого микроскопа (СЗМ). Для этого в уже имеющиеся аппаратные схемы, предложенные в 
работе 1, нужно добавить фазосдвигающие элементы, задающие фазовые сдвиги Да, 
Ас (или Аа, Ао). 

Вместо аппаратного возведения в степени т и п синусоидального сигнала с помощью 
умножителей, как это было предложено в работе 1, можно согласно (3) просто суммиро- 
вать синусоидальные и косинусоидальные сигналы с кратными частотами. Например, для 
компенсации искажений, вызываемых петлей гистерезиса Лист (т=3, п=1), показанной на 
Рис. 37а, потребуется два косинусоидальных сигнала с частотами фи За, элемент, сдви- 
гающий фазу на четверть периода, три усилителя с коэффициентами усиления З/4а, а/4, 
Ь; и два суммирующих усилителя. 

Согласно (7) для компенсации искажений петли Лист потребуется два косинусоидаль- 


ных сигнала с частотами @ и З@ элемент, сдвигающий фазу на величину 


ф.=агсап(46/За), два усилителя с коэффициентами усиления Апи=_/9а?/16 + 6? , Ать=а/4 


и суммирующий усилитель.? 


У. РЕЗЮМЕ 

Выполнена доработка существующей модели петли гистерезиса, построенной на па- 
раметрических уравнениях; уточнено несколько формул, полученных ранее. В результате 
доработки в несколько раз увеличена точность аппроксимации. Показано, что порождаю- 
щую функцию петли гистерезиса можно представить в виде суммы нерасщепленной пет- 
ли и кривой расщепления, в виде частотного спектра и в экспоненциальном виде. 

Выведена общая формула, с помощью которой можно строить кусочно-линейные пет- 
ли гистерезиса типа Люфт и Неидеальное Реле, а также их многочисленные разновидно- 
сти, широко используемые при построении упрощенных моделей гистерезисных явлений. 
Продемонстрирована возможность образования/декомпозиции различных петель гисте- 
резиса. Для ряда типов петель результаты образования/декомпозиции дают основание 
предположить одновременное существование в рассматриваемой системе не одного, а 
двух различных гистерезисных процессов, наложенных друг на друга. Получено несколь- 
ко новых формул, описывающих различные виды одинарных, двойных и тройных петель. 

Помимо представленных параметрических уравнений параллелограмма и прямо- 
угольника, найдены параметрические уравнения ромба, квадрата, правильного шести- 
угольника и правильного восьмиугольника.° Предложено более общее выражение для 
вычисления площади петли гистерезиса, служащей для оценки потерь энергии в пьезо- 
электрике/ферромагнетике. В ходе исследования обнаружено несколько тождеств, свя- 
зывающих биномиальные коэффициенты.? 


Разработанная модель петли гистерезиса особенно удобна при решении задач ими- 
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тационного моделирования циклически работающих приборов, включающих в себя зве- 
нья с гистерезисом. Кроме того, модель позволяет простыми аппаратными средствами 
сформировать выходные сигналы х(й, у(й, обеспечивая таким образом легкую аппарат- 


ную реализацию как прямых, так и обратных петель гистерезиса. 


ДОПОЛНИТЕЛЬНЫЙ МАТЕРИАЛ 

Дополнительный материал включает 27'р-архив рабочих листов Маткада 20011, в кото- 
рых детально рассмотрены все аспекты исходной и улучшенной параметрических моде- 
лей петли гистерезиса. Те, у кого нет программы Маткад, могут воспользоваться 
читаемыми рабочими листами Маткада в виде РОЕ-документа, который прилагается. Из- 
за ограничения на размер статьи в ней представлены только наиболее часто встречаю- 
щиеся на практике петли гистерезиса. Если требуемая петля отсутствует в статье, то её 


поиск стоит продолжить в дополнительном материале. 


Работа выполнена при финансовой поддержке “Фонда содействия развитию малых 
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выражает благодарность О. Э. Ляпину, О. \\/. \/аааей, о. В. Объедкову и Л. Б. Шаровой за 
критическое прочтение рукописи и проверку дополнительных материалов, доц. 


Е. А. Фетисову за поддержку и стимулирование работы. 
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